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I.  INTRODUCTION 


The  in-flight  characteristics  of  an  artillery  shell  are  of  major 
importance  to  the  shell  designer,  ballistician  and  ultimately  the  artillery 
field  commander  whose  mission  is  deployment  of  timely  and  accurate  fire 
power.  The  aerodynamic  properties  of  artillery  shell,  such  as  pitching  moment 
and  Magnus  moment  are  critical  to  the  stability  of  shell  which  in  turn  signi¬ 
ficantly  affect  accuracy  and  time  of  flight.  The  capability  to  determine  the 
aerodynamics  of  shell  is  required  over  a  wide  range  of  flight  regimes  since, 
depending  on  initial  launch  velocities,  artillery  shell  are  subject  to 
subsonic,  transonic,  and  supersonic  flight.  Projectile  aerodynamics  over 
these  various  flight  regimes  have  been  found,  in  some  cases,  to  change  by  an 
order  of  magnitude.  The  solution  techniques  utilized  must  therefore  be 
capable  of  computing  these  changes. 

A  concentrated  theoretical  and  experimental  research  program  has  been 
ongoing  at  8RL  in  order  to  develop  the  predictive  capabilities  required  for 
determining  projectile  aerodynamics.  Supersonic  computations  using  combined 
inviscid  flow  field  and  boundary  layer  techniques  have  been  developed  by 
Sturek1,  et  al.,  for  cone-cylinder  and  ogive-cylinder  configurations.  Recent 
results  have  been  obtained  in  supersonic  flow  over  a  typical  boattailed  pro¬ 
jectile  by  Schiff  and  Sturek2  using  modern  computational  techniques  for 
solving  the  thin-layer  Navier-Stokes  equations. 

Inviscid  transonic  computational  results  have  been  obtained  by 
Reklis3,  et  al.,  for  a  secant-ogive-cylinder-boattail  shape.  The  inviscid 
techniques  give  fair  results  for  pitch  plane  aerodynamic  coefficients  at  small 
angle  of  attack.  The  ability  to  compute  the  Magnus  effect,  however,  relies  on 
an  accurate  computation  of  the  viscous  boundary  layer.  Techniques  which  have 
been  applied  in  supersonic  flow  for  combining  inviscid  and  boundary  layer 
methods  have  not  been  fully  established  for  transonic  flow.  These  methods, 
which  have  shown  good  results  for  ogive-cylinders  in  supersonic  flow  at  low 
angle  of  attack,  are  not  accurate,  however,  in  modeling  the  severe  flow 
expansion  in  the  vicinity  of  surface  discontinuities  such  as  those  that  occur 
at  the  cylinder-boattai 1  junction. 


2.  Sturek,  W.  B.,  et  al.,  "Computations  of  Magnus  Effects  for  a  Yawed, 

Spinning  Body  of  Revolution AIAA  Journal,  Vol.  16,  No.  7,  July  1978,  pp. 
687-692. 

2.  Schiff,  L.  B.,  and  Sturek,  W.  B.,  " Numerical  Simulation  of  Steady 
Supersonic  Flew  Over  Cone  Ogive-Cylinder-Boattail  Body,"  AIAA  Paper 
No.  80-0066,  14-16  January  1980. 

3.  Reklis,  R.  P.,  Sturek,  W.  B. ,  and  Bailey,  F.  L. ,  "Computation  of  Transonic 
Flow  Past  Projectiles  at  Angle  of  Attack U.S.  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  Technical  Report  ARBRL-TR-02139,  February  1979. 

AD  A069106. 
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The  solution  of  the  thin- layer  Navier-Stokes  equations,  which  allows  for 
the  simultaneous  computation  of  the  invlscid  and  viscous  regions,  eliminates 
the  need  for  matching  two  different  solutions.  Additionally,  since  all  three 
momentum  equations  are  retained,  the  ability  to  compute  in  regions  of 
separated  flow  is  achieved.  This  paper  describes  the  governing  three- 
dimensional  thin-layer  Navier-Stokes  equations  used  for  computing  flow  over 
projectile  shapes  at  angle  of  attack.  Secondly,  the  generalized  ax i symmetric 
formulation  used  for  computations  at  a  =  0°  will  be  described.  A  description 
of  the  numerical  algorithm  and  results  will  follow.  Computational  results  and 
limited  experimental  data  will  be  presented  for  a  secant-ogive-cylinder- 
boattail  projectile  shape  (SOCBT)  at  a  =  0°.  Additional  computations  have 
been  obtained  for  the  same  projectile  shape  at  a  =  2°.  Computational  results 
will  also  be  shown  for  a  ring  airfoil  shape  thus  demonstrating  the  general 
geometry  capability  of  the  present  numerical  scheme. 

II.  GOVERNING  EQUATIONS 


The  general  three-dimensional  thin-layer  Navier-Stokes  equations,  used 
for  all  cases  where  a  >  0°,  are  described  in  Section  IIA.  The  thin-layer 
generalized  axisymnetric  equations,  which  are  a  special  case  in  the  3-0 
equations,  are  described  in  Section  I  IB. 


Three-Dimensional  Equations 


The  transformed  three-dimensional  thin-layer  Navier-Stokes  equations  in 
non-dimensional  and  strong  conservation  law  form  are  written  as4 


9xq  +  3^E  +  +  3^G  =  Re  3^S 


where  general  coordinate 

S  =  S(x,y,z,t)  -  longitudinal  coordinate 
n  =  n(x,y,z,t)  -  circumferential  coordinate 
c,  =  c(x,.y,z,t)  -  near  Normal  coordinate 
t  =  t  -  time 


4.  Pulliam,  T.  H.,  and  Steger,  J.  L.,  "On  Implicit  Finite-Difference 
Simulations  of  Three-Dimensional  Flow,"  AIAA  Journal,  Vol .  18,  No.  2, 
February  1980,  pp.  159-167. 


are  used  and 
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The  velocities  in  the  £,n,£  coordinates  are 

U  =  £f  +  £u  +  £  v  +  £  w 

v  A  J  Z 

v  =  nt  +  nxii  +  nyv  +  ry*  (2) 

W  =  £t  +  £  u  +  £  v  +  £  w 

represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  components  (u,v,w)  are  retained  as  the  dependent 
variables  and  are  nondimensionali zed  with  respect  to  a^  (the  free  stream  speed 

of  sound).  The  local  pressure  is  determined  using  the  relation 

p  =  (y  -  l)[e  -  .5p(u2  +  v2  +  w2)]  (3) 

where  y  is  the  ratio  of  specific  heats,  density  (p),  is  referenced  to  px  and 

total  energy  (e)  to  pM  aM2.  The  additional  parameters  are  (<)  the  coefficient 

of  thermal  conductivity,  (w)  the  dynamic  viscosity,  (Re)  the  Reynolds  number, 
(Pr)  the  Prandtl  number,  and  (x)  which  through  the  Stokes  hypothesis 
is  (-2/3) u. 

The  metric  terms  of  Equation  (1)  are  defined  from 
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The  “thin-layer"  approximation4"7  used  here  requires  that  all  body  sur¬ 
faces  be  mapped  onto  ?  =  constant  planes  and  Re  »  1.  Essentially,  all  the 
viscous  (second  derivative)  terms  in  the  coordinate  directions  (here  taken  as 
£  and  n)  along  the  body  surface  are  neglected  while  terms  in  the  e  or  the  near 
normal  direction  to  the  body  are  retained.  This  approximation  is  used 
because,  due  to  computer  speed  and  storage  limitations,  fine  grid  spacing  can 
only  be  provided  in  one  coordinate  direction  (usually  taken  as  the  near  normal 
direction)  and  the  grid  spacing  available  in  the  other  two  directions  is 

usually  too  coarse  to  resolve  the  viscous  terms.  For  the  type  of  problems 
currently  under  investigation,  i.e.,  projectiles  at  low  angles  of  attack,  with 
no  strong  cross-flow  separation,  these  approximations  are  considered  valid. 

B.  Generalized  Axisymmetric  Equations. 

The  thin-layer  general i zed- axisymmetric  equations  are  obtained  from  the 
three-dimensional  equations  by  making  use  of  two  restrictions:  (1)  all  body 
geometries  are  of  an  axisymmetric  type;  (2)  the  state  variables  and  the 

contravariant  velocities  do  not  vary  in  the  circumferential  direction.  In 

what  follows,  the  a^F  term  of  Eq.  (1)  will  be  reduced  to  the  source  term  of 

the  generalized  axisymmetric  equations. 

A  sketch  of  a  typical  axisymmetric  body  is  shown  in  Figure  la.  In  order 
to  determine  the  ci rcumferential  variation  of  typical  flow  and  geometric 
parameters,  we  first  establish  correspondence  between  the  inertial  Cartesian 
coordinates  (x,y,z)  (to  which  the  dependent  variables  are  referenced),  the 
natural  inertial  cylindrical  coordinates  (x,<p,R),and  the  transformed  variables 
U,n,t).  The  choice  of  the  independent  variables  K,  n,  c  is  restricted,  as 

shown  in  Figure  lc,  insofar  as  n  must  vary  as  <t> ,  i.e.,  $  =  Cn  (where  C  is  a 

constant).  From  the  views  shown  in  Figure  1,  the  relationship  between  the 
coordinate  systems  are  observed  to  be 


5.  Nietubicz,  C.  J.,  Pulliam ,  T.  H.,  and  Steger,  J.  L.,  "Numerical  Solution 
of  the  Azimuthal-Invariant  Thin-Layer  Navier-Stokes  Equations U.S.  Army 
Ballistic  Research  Laboratory,  ARRADCOM,  Te~hniaal  Report  ARBRL-TR-02227, 
March  1980.  AD  AO 85716. 

6.  Steger,  J.  L.,  "Implicit  Finite  Difference  Simulation  of  Flow  About 
Arbitrary  Geometries  with  Application  to  Airfoils,"  AIAA  Paper  No.  77-665, 
June  1977. 

7.  Baldwin,  B.  S.,  and  Lomax,  H.,  "Thin  Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flows,"  AIAA  Paper  No.  78-257,  January  1978. 
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<t>  =  Cn 
x  =  x(^,4,x) 
y  =  R(^,<.,x)sin$ 
z  =  R(S,c,t)cos$ 

where  <$>  =  and  the  Cartesian  and  cylindrical  coordinates  are  related  in 
the  usual  way.  Note  that  x  and  R  are  general  functions  of  only  £,  c,  and  t. 

Evaluating  the  metric  terms  given  the  above  assumptions  and  substituting 
,  in  Equation  (1)  the  resulting  unsteady  thin-layer  generalized  axisymmetric 
equation5  can  be  written  as 

3rq  +  3^E  +  3^6  +  H  =  Re'^S  (6) 

where 

0 
0 

h  =  j-1^  Pv[Rc(u-ct)  +  Rt(w-ct 
-pVR»n(V-nt)  -  p/(Rvn) 

0 

is  the  resultant  source  term  which  has  replaced  3r)E  of  Equation  (1). 

Equation  (6)  contains  only  two  spatial  derivatives  but  does  retain  all 
three  momentum  equations  thus  allowing  a  degree  cf  generality  over  the 
standard  axisymmetric  equations.  In  particular,  the  circumferential  velocity 
is  not  assumed  to  be  zero  allowing  then  computations  for  spinning  projectiles 
or  swirl  flow  to  be  accomplished. 

The  numerical  algorithm  used  for  both  Equations  (1)  and  (6)  is  a  fully 
implicit,  approximately  factored  finite  difference  scheme  as  analyzed  by  Beam 
and  Warming8.  Additional  details  of  the  numerical  method,  algorithm,  and 
boundary  conditions  for  each  formulation  can  be  found  in  References  4  and  5, 
respecti  vel.y. 


8.  Beam ,  R.,  and  Warming,  R.  F.,  "An  Implicit  Factored  Scheme  for  the 

Compressible  Navier-Stokes  Equations ,"  AIAA  Paper  No.  77-645,  June  1977. 


1. 


III.  RESULTS 


A  series  of  computations  has  been  obtained  for  a  3  caliber  secant-ogive 
nose,  2  caliber  cylinder,  and  1  caliber,  7°  boattail  shape.  The  generalized 
axi symmetric  formulation  was  used  for  all  cases  where  «  =  0°  and  the  3-D 
Navier-Stokes  Pulliam-Steger  code  was  used  for  u>  0°.  The  results  are 
presented  in  the  form  of  contour  plots,  surface  pressure  distribution,  and 
velocity  profiles.  Computations  have  been  performed  for  a  ring  airfoil  shape, 
the  results  of  which  show  an  interesting  shock  pattern  being  developed  as  a 
function  of  Mach  number. 

A.  Secant-Ogi  ve-Cylinder-Boattail ,  a=  0°. 

A  computational  finite  difference  mesh  of  78  longitudinal  points  by  50 
normal  points  was  used  for  these  calculations.  Grid  points  were  clustered  in 
the  vicinity  of  the  body  surface  discontinuities,  which  exist  at  the  ogive- 
cylinder  and  cylinder-boattail  junction.  The  grid  points  in  the  normal 
direction  were  exponentially  stretched  away  from  the  surface  with  approxi¬ 
mately  25  points  located  within  the  boundary  layer.  The  Reynolds  number  for 
all  cases  was  0.72  x  106  based  on  body  diameter  and  free  stream  velocity.  The 
initial  computation  was  for  a  free  stream  condition  of  Mach  =  0.8.  Additional 
runs  were  made  by  incrementing  the  free  stream  Mach  number  and  using  the 
previous  converged  solution  as  a  starting  condition.  The  results  are  first 
shown  in  the  form  of  Mach  contours  in  Figure  2.  The  regions  of  subsonic  and 
supersonic  flow  are  identified  and  the  coalescence  of  the  Mach  lines  to  the 
right  of  the  supersonic  region  represents  the  position  of  the  shock. 
Initially  at  M^  =  0.8  (Figure  2a),  only  a  small  disturbance  is  felt  in  the 

free  stream.  As  the  Mach  number  is  increased  to  0.92  (Figure  2b),  shocks  are 
seen  to  occur  on  the  cylinder  and  boattail.  With  a  further  increase  in  Mach 
number,  the  regions  of  supersonic  flow  have  grown  and  the  shocks  have  moved 
well  downstream.  A  further  increase  in  the  free  stream  condition  to  Mm  =  1.1 

(Figure  2f)  shows  the  beginning  of  a  bow  shock  at  the  nose  and  expansion  fans 
are  evident  at  the  surface  discontinuities. 

A  more  critical  look  at  the  computational  results  is  shown  in  Figure  3 
where  the  surface  pressure  coefficient,  Cp,  is  shown  as  a  function  of  axial 

position.  The  computational  results  are  indicated  by  a  solid  line  and  the 
circles  are  experimental  data  obtained  by  Kayser9.  Overall  the  comparisons 
show  generally  good  agreement,  however,  some  discrepancy  is  apparent  near  the 
expansions  (decrease  in  Cp)  and  shocks  (sharp  increase  in  Cp)  which  indicates 

the  need  for  additional  mesh  points  in  these  areas.  Additionally,  the 
changing  nature  of  the  pressure  distributions  indicate  that  the  varying  shock 


.9.  Kayser,  l.  D.,  and  Whiton,  F.,  "Surface  Pressure  Measurements  on  a 

Boattailed  Projectile  Shape  at  Transonic  Speeds BRL-MR  to  he  published. 


structure  has  been  predicted  quite  well.  An  accurate  computation  of  the 
pressure  distribution  is  important  since  the  integration  of  surface  pressures 
is  a  major  contributor  in  determining  the  aerodynamic  coefficients.  This  is 
especially  apparent  in  the  transonic  flight  regime  where  the  magnitude  of  the 
coefficients  can  change  by  as  much  as  100%.  For  the  series  of  computations 
obtained  from  =  0.8  through  =  1.1,  the  surface  pressures  were  integrated 

to  determine  the  aerodynamic  wave  drag.  The  results  are  shown  in  Figure  4 
together  with  the  available  experimental  data.  Excellent  agreement  is  shown 
for  both  the  drag  rise  and  magnitude  of  the  wave  drag  in  the  critical  Mach 
number  regime. 

The  dynamic  stability  of  shell  is  one  area  of  concern  when  designing  new 
shell  or  modifying  existing  ones.  The  Magnus  moment,  which  affects  the 
dynamic  stability,  is  a  viscous  phenomenon  occurring  at  a  >  0°.  Therefore  an 
accurate  representation  of  the  viscous  portion  of  the  flow  field  is  crucial  to 
computing  the  Magnus  moment.  As  an  initial  examination  of  the  boundary  layer 
in  transonic  flow,  with  its  associated  shock  interaction,  computations  were 
performed  at  =  0.97,  Reg  =  2.6  *  106  and  a  =  0°  for  a  similar  projectile 

shape.  A  comparison  of  the  computational  velocity  profiles  with  experimental 
data  obtained  by  Danberg10  is  shown  in  Figure  5.  The  results  are  shown  for 
three  axial  stations  in  the  vicinity  of  the  boattail  (X/D  =  5.3).  The  pro¬ 
files  compare  very  well  at  X/D  =  5.05  and  5.61,  which  are  stations  approxi¬ 
mately  .3  in  front  of  and  behind  the  boattail  corner.  The  profile  at  X/D  = 
5.36,  which  is  in  an  area  of  severe  expansion,  is  also  considered  quite 
good.  The  adequate  resolution  of  the  boundary  layer  which  in  this  case 
includes  a  region  of  shock  boundary  layer  interaction  is  very  encouraging  for 
future  work  in  computing  the  Magnus  effect. 

B.  Secant-Qgive-Cyl i nder-Boattai 1 ,  a  =  2°. 

A  new  finite  difference  mesh  consisting  of  60  longitudinal  points,  28 
normal  points,  and  20  points  in  the  circumferential  direction  was  designed  for 
computations  at  angle  of  attack.  Clustering  of  the  longitudinal  points  was 
maintained  in  the  vicinity  of  the  expansion  similar  to  the  a  =  0°  cases.  The 
computed  leeward  and  windward  pressure  distributions  are  shown  in  Figure  6. 
The  pressure  on  the  ogive-nose  is  seen  to  be  higher  on  the  windward  side,  thus 
generating  a  positive  or  upward  force  on  the  projectile.  The  boattail  region, 
however  shows  a  reversal  with  the  higher  pressure  now  occurring  on  the  leeward 
ray  thus  causing  a  negative  or  downward  force.  This  resulting  couple  about 
the  center  of  gravity  is  the  main  contributor  to  the  critical  aerodynamic 
behavior  in  the  transonic  regime.  This  effect  can  additionally  be  seen  in  the 
Mach  contours  of  Figure  7.  The  asymmetry  in  the  flow  field  is  apparent 


10.  Danberg,  J .,  Wind  tunnel  data,  to  be  published  as  a  University  of 
Delaware  Report  (Private  communications) . 
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between  the  wind  and  lee  side  and  more  importantly  the  boattail  shock  is  seen 
to  be  further  aft  on  the  wind  side. 

Velocity  distributions  obtained  from  the  3-0  computations  are  shown  in 
Figure  8  for  three  axial  stations.  The  roll  angle  is  located  at  the  top  of 
each  profile  with  <j>  =  0°  and  180°  corresponding  to  the  lee  and  wind  side, 
respectively.  There  are  no  experimental  data  for  this  case,  however,  the 
profiles  do  show  a  smooth  variation  with  roll  angle  and  become  more  full  as 
the  lee  side  is  approached.  Additionally,  looking  at  the  <p  =  0°  plane  for  the 
three  axial  stations,  the  velocity  is  seen  to  increase  near  the  boattail 
corner  (X/D  =  5.0)  where  the  flow  is  expanding  and  then  decrease  on  boattail 
(X/D  =  5.6). 

C.  Hollow  Projectile,  m  =  0° . 

Of  current  interest  in  shell  design  is  the  utilization  of  hollow  projec¬ 
tiles  which  have  the  characteristic  of  "flat"  trajectories.  A  shape  of  this 
type,  known  as  the  ring  airfoil,  has  been  type  classified  and  is  currently 
used  as  an  anti-riot  device.  In  order  to  demonstrate  the  capability  to 
compute  flows  about  hollow  projectile  shapes  of  current  Army  interest,  the 
generalized  axi symmetric  Navier-Stokes  code  was  used  for  computations  of  a 
ring  airfoil  shape  at  a  =  0°.  A  cross  section  of  the  actual  shape  is  shown  at 
the  bottom  of  the  Cp  plot  of  Figure  9a.  Inviscid  results  are  presented  in 

Figures  9a,  b,  c,  and  d  for  Mw  =  0.4,  0.7,  0.8,  and  0.9.  In  all  cases  the 

pressure  distribution  is  plotted  for  the  internal  and  external  surfaces  using 
a  solid  line  and  dashed  line,  respectively. 

The  critical  pressure  coefficient,  Cp*.  is  the  value  of  the  pressure 

coefficient  at  sonic  velocity.  The  flow  velocity  in  regions  with  pressure 
greater  than  Cp*  is  subsonic  and  in  regions  where  the  pressure  is  less  than 

Cp*  the  flow  is  supersonic.  For  =  0.4  (Figure  9a),  all  values  of  Cp  are 

greater  than  Cp*  indicating  the  flow  over  both  the  internal  and  external 

surfaces  is  subsonic.  However,  as  the  Mach  number  is  increased  to  0.7,  the 
upper  surface  is  shown  to  develop  supersonic  flow  and  a  shock  wave  while  the 
internal  flow  remains  subsonic.  Increasing  the  Mach  number  still  further 
to  M^  =  0.8,  the  shock  waves  are  now  seen  to  exist  on  both  the  external  and 

internal  surfaces.  A  final  solution  at  M^  =  0.9  shows  that  the  shock  wave  has 

moved  to  the  trailing  edge  indicating  supersonic  flow  over  most  of  the 
internal  and  external  surfaces. 


IV.  SUMMARY 

An  implicit,  finite  difference  technique  has  been  used  to  solve  the  thin- 
layer  Navier-Stokes  equations.  The  three  dimensional  and  generalized  axisym- 
metric  equations  have  been  presented.  Solutions  are  obtained  for  a  secant 
ogive  boattail  and  hollow  projectile  shape. 
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Comparisons  of  experimental  and  computed  pressure  distributions  were 
obtained  from  Mach  =  0.94  to  1.1.  The  comparisons  indicated  excellent  agree¬ 
ment  throughout  the  entire  Mach  number  range.  The  movement  of  the  shock 
position  corresponding  to  changes  in  the  free  stream  condition  can  be  clearly 
seen  from  the  Mach  number  contour  plots. 

The  Navier-Stokes  computations,  which  implicitly  model  the  viscous  and 
inviscid  flow  field,  indicate  excellent  agreement  with  experimental  velocity 
profile  data  for  the  a  =  0  case.  This  result  is  very  encouraging  since 
accurate  prediction  of  the  turbulent  boundary  layer  development  is  critical  to 
the  determination  of  the  Magnus  force. 

Results  have  been  obtained  for  a  projectile  shape  at  «  =  2.0,  Mach  = 
0.96.  The  general  features  of  the  flow  field  have  been  obtained.  The  Mach 
number  contour  plot  shows  the  asymmetric  shock  pattern  which  exists  on  the 
boattail  and  pressure  plots  of  wind  versus  lee  side  indicate  the  correct 
trends.  A  comparison  of  these  results  to  experimental  surface  pressure  and 
velocity  profile  data  are  required,  however,  to  assess  the  accuracy  of  the 
three  dimensional  computations. 

The  general  capability  to  compute  internal  flows  has  been  demonstrated  by 
solving  the  inviscid  flow  field  about  a  ring  airfoil  projectile.  Computations 
were  obtained  for  Mach  =  0.4  to  0.9.  A  significant  change  in  the  flow  field 
has  been  found  to  occur  throughout  this  regime.  The  Mach  contours  presented 
indicate  reqions  of  both  subsonic  flow  and  mixed  subsonic/supersonic  flow. 

The  results  presented  here  provide  a  significant  measure  of  confidence  in 
using  these  new  computational  techniques  for  calculating  projectile  flow 
fields. 
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PROJECTILE -TYPE  »O0Y  X  »  CONST.  PLANE 
(a)  (b) 


♦  «Ci» 

X  »  X(f,£,r) 

Y  -R(e,Ctr)»in* 
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Figure  1.  Axisymmetric  body  and  coordinate  system 
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Figure  Zb.  Mach  contours  for  SOCBT,  a 
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Figure  2f.  Mach  contours  for  SOCBT,  a=0,  M= 1.1 
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Figure  4.  Wave  drag  for  SOCBT,  a  =  0.  Theory  and 
experiment. 
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Figure  5.  Comparison  of  theoretical  and  experimental 
velocity  profiles  for  SOCBT. 
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Figure  9a.  Inviscid  transonic  flow  over  a  ring  airfoil 
projectile,  M  =  0.4,  a  =  0° 
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Figure  9b.  Inviscid  transonic  flow  over  a  ring  airfoil 
projectile,  M  =  0.7,  a  =  0° 
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Figure  9c.  Inviscid  transonic  flow  over  a  ring  airfoil 
projectile,  M  =  0.8.  a  =  0° 


Figure  9d.  Inviscid  transonic  flow  over  a  ring  airfoil 
projectile,  M  =  0.9,  a  =  0° 
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LIST  OF  SYMBOLS 


a 

C 

Co 

CP 

Cr 

D 

e 

E  ,F  ,G  ,H  ,S  ,q 
J 
M 
P 

Pr 

R 

Ren 


Re 

t 


u,v,w 

u.v.w 


x,y  a 


r 

< 

M 


-  speed  of  sound 

-  proportional  constant  (see  Equation  5) 

2 

-  pressure  coefficient,  (p-P J/ lh 

-  specific  heat  (at  constant  pressure) 

-  cord  lenqth  (see  Fiqure  R) 

-  body  diameter  (57.15mm) 

-  total  energy  per  unit  volume/p^a^ 

-  flux  vector  of  transformed  Navier-Stokes  equations 

-  transformation  Jacobian 

-  Mach  number 

2 

-  pressure/p^u^ 

-  Prandtl  number,  m  c  /< 

9  CO  p  '  CD 

-  body  radius 

-  Reynolds  number 

-  Reynolds  number,  Rep/M^ 

-  physical  time 

-  Cartesian  velocity  components/a^, 

-  Contravariant  velocity  component s/a^ 

-  physical  Cartesian  coordinates 

-  anqle  of  attack 

-  ratio  of  specific  heats 

-  coefficient  of  thermal  conductivity 

-  coefficient  of  viscosity 
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Superscript: 

★ 


Subscript: 


-  transformation  coordinates  in  axial,  circumferential 
radial  directions 

-  density/p^ 

-  transformed  time 

-  circumferential  angle 

-  critical  value 

-  free  stream  conditions 
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